Relations de récurrence linéaires, primitivité et loi de 

Benford 

Hugues Deligny et Paul Jolissaint 
18 août 2010 
O- 
(Ni ; 1 Introduction 

Nous nous intéressons ici à la loi de Benford pour des suites (a„) n >i C IR qui satisfont 
"^ ! une relation de récurrence d'ordre k de la forme : 

(*) O-n+k = Ck-lO-n+k-1 + Ck-20 n +k~2 + ■ ■ ■ + ClO n+ l + Coa n 

avec les conditions initiales Oj > pour tout 1 < i < k. On supposera la plupart du 

temps que les coefficients q sont positifs ou nuls, et que cq > 0, de sorte que la relation de 

récurrence soit effectivement d'ordre k. On va présenter des conditions suffisantes sur les 

Ci pour qu'il existe (au moins) un entier b > 2 tel que la suite (a n ) n >i satisfasse la loi de 

Benford en base b. 

Rappelons qu'étant donné un tel entier b > 2, une suite (w n ) n >i C M + \ {0} satisfait la 

(D ■ loi de Benford en base b si, pour tout t 6 [1, b), on a : 
> ■ 

S? 1 r |{l<n<iV : M h (u n )<t}\ 

S ' i im AT = log ^)' 

in ; 

où Mb(u n ) désigne la mantisse de u n , c'est-à-dire l'unique élément de [1, b) tel que u n = 
Mb(u n ) ■ b m avec m G Z (cf [5], définition 2.1). La définition ci-dessus généralise le cas 
classique qui correspond à b = 10 (dix) et qui affirme que le premier chiffre significatif 
di(u n ) de u n satisfait la loi de probabilité (cf remarque (2) de [S]) : 
'O ' pour tout chiffre 1 < d < 9, on a 

\{l<n<N:d 1 (u n ) = d}\ ( 1 

lim = log 10 1 + 3 

D'après les théorèmes principaux de [1] et [5], une suite (u n ) C M + \ {0} satisfait la loi 
de Benford en base b dès qu'elle remplit les deux conditions suivantes : 

u 

(a) il existe des nombres réels a>0,p>0et/i tels que lim — — = a ; 

n— >oo fit 1 n n 

(b) log fe (p) est irrationnel. 

Si c'est le cas, si Q(x) est un polynôme non constant à coefficients entiers et Q{x) > 1 
pour x assez grand, la sous-suite {uQ^ n )) n >i satisfait également la loi de Benford en base b. 

L'exemple suivant découle immédiatement des deux conditions ci-dessus, et il ne semble 
pas avoir été découvert auparavant : 

1 



O 
O 



Exemple 1 Soit (j) n )n>i la suite croissante des nombres premiers. Comme cela a été ob- 
servé dans [3j, la suite (p n ) ne satisfait pas la loi de Benford. En revanche, choisissons deux 
entiers £ > 2 et b > 2 tels que log 6 (£) soit irrationnel. Alors la suite (pp) satisfait la loi de 
Benford en base b. En effet, par le théorème des nombres premiers, on a 

lim -, r = 1. 

n-i-oo nln(n) 

En passant à la sous-suite (£ n ) n >i C N*, on a 

lim — — = 1 et donc lim — — = \n(£). 

™->oo n£ n m(£) n->oo n£ n 

Dans le cas d'une suite (o n ) n >i qui satisfait une relation de récurrence du type (*), 
(a n ) n >i s'exprime à l'aide des racines du polynôme caractéristique de la relation de récur- 
rence, c'est-à-dire le polynôme 

p(x) = X k — Ck-\X k ~ l — ... — C\X — c . 

Plus précisément, écrivons p(x) = (x — Ci) Ml • • • (x — Cm) Mm où Ci> • • • ; Cm G C sont les 
zéros distincts de p(x) de multiplicités respectives Lii, . . . , \x m > 1. Alors (o n ) n >i est une 
combinaison linéaire des k suites (?/C n )«>i pour < £ < Llj et 1 < j < m : il existe des 
constantes a,j/ dépendant des conditions initiales oi, . . . , Ofe telles que 

3=1 1=0 

Dès lors, si p(x) admet une racine positive simple p telle que \Ç\ < p pour toute autre 
racine ( de p(x) et si le coefficient de p n est positif dans la décomposition ci-dessus, alors 
la suite (a n ) n >i satisfait la condition (a). 

Il nous a semblé intéressant de donner des conditions suffisantes sur les coefficients de 
la relation de récurrence (*) plutôt que sur les racines du polynôme caractéristique de 
la relation de récurrence. Nous mettrons plus particulièrement l'accent sur le cas où les 
coefficients q sont positifs ou nuls. 

Le premier résultat de notre étude est : 
Théorème 2 Soient cq, Ci, . . . , Ck-\ des nombres réels. On pose 

I = {l<j<k-l: Cj ^0}, 
et on suppose que 1^0. Enfin, on associe à la suite Cq, . . . , c^-i le polynôme 

p(x) = X k — Ck-\X k ~ l — ... — C\X — Cq. 

(1) Si les Ci sont positifs ou nuls et si Cç, > 0, le polynôme p admet une unique racine 
positive p et \(\ < p pour tout autre zéro C, de p (compté avec multiplicité) si, et 
seulement si pgcd(J U {k}) = 1. 

(2) Si Cfc_i > 2 et si Ck-i > Y2jZq \ c j\ + 1? alors le polynôme p admet une racine simple 
Ck-\ — 1 < p < Cfe^i + 1, et toute autre racine ( de p satisfait \(\ < 1. 



(3) Supposons que les Ci satisfassent l'une des deux conditions ci-dessus, soit b > 2 un 
entier tel que log b (p) soit irrationnel, et soit une suite (a„) n >i C M + \{0} qui satisfait 
la relation de récurrence 

dn+k — Ck-ia n+ k-i + . . . + C\a n+ i + c^a n . 

Si les Ci satisfont la condition (1) ou la condition (2), et dans ce dernier cas si de 
plus (o n ) ne tend pas vers 0, alors la suite (a n ) n >i satisfait la loi de Benford en 
base b. Enfin, il en est de même de toute sous-suite de la forme (a<2( n )) où Q(x) 
désigne un polynôme non constant à coefficients entiers tel que Q(x) > 1 pour tout 
x suffisamment grand. 

Remarque. La condition supplémentaire a n -/ï dans la troisième affirmation du théorème 

ci-dessus est indispensable. En effet, la suite a n = 10 _n tend vers 0, satisfait la relation de 

récurrence 

31 3 

a n+2 — Tpf a n+1 — ï~n an 

dont le polynôme caractéristique est p(x) = (x — 3)(x — 1/10) ; ses coefficients satisfont la 
condition (2) du théorème, mais la suite ne satisfait pas la loi de Benford pour b = 10, bien 
que log 10 (3) soit irrationnel. 

On verra que lorsque les q sont positifs ou nuls, le polynôme caractéristique p(x) = 
x h — Ck-iX k ~ 1 — ... — C\X — Co admet toujours une racine dominante p > 0, c'est-à-dire 
telle que \(\ < p pour toute autre racine ( ; regardons maintenant le cas où il admet h > 1 
racines (i = p, . . . , (h telles que \Q\ = p. On note encore I = {1 < j < k : Cj > 0}, et pour 
tout m G {0, . . . , h — 1}, on pose 

Im — {j £ I '■ j — ^(mod h)} 

de sorte que (I m )o< m <h constitue une partition de /. 

Théorème 3 Soient cq > 0, c\ > 0, . . . , Ck-i > 0, et h comme ci-dessus. Alors h = 
pgcd(J U {k}). De plus, toute suite (a n ) n >i qui satisfait la relation de récurrence 

O-n+k — Ck-lO-n+k-l + Cfc-2Ûn+fe-2 + • • • + Cl^n+l + CoGt„ 

avec des conditions initiales positives ou nulles ai,...,dk est réunion des h sous-suites 
((dm+hn)n>i)o<m<h- Enfin, pour tout < m < h, on a l'alternative suivante : 

- Si {j G I m : dj > 0} 7^ 0, alors linij^oo am + n hn existe et est positive, et la suite 
(d m+ hn)n>i satisfait la loi de Benford en base b pour tout b tel que log b (p) ^ Q. 

- Si {j G I m '■ dj > 0} = ; alors a m+ h n = pour tout n > 1. 

Exemple 4 Le théorème 2 s'applique au cas des suites de Fibonacci d'ordre k > 2 : on 
choisit des nombres réels ai, . . . , a& > arbitraires, et on définit (a n ) n >i par 

d n +k — On+fc-1 + d n+ k-2 + • • • + d n+ \ + d n . 

Notons que la racine dominante p > est irrationnelle grâce à l'observation suivante : 

Proposition 5 Soit k > 2 un entier et soient a, . . . , Ck~i G N, et p(x) := x k — Ck-iX — 
... — C\X — 1. Si p > est une racine rationnelle de p(x), alors p — 1. 



Preuve. Écrivons p = - avec p, q G N* et pgcd(p, q) — 1. Alors, en utilisant le fait que p est 
une racine de p(x), on obtient 

P k = Cfe-i/ _1 g + c k ^ 2 y k - 2 q 2 + ... + c x pq k ~ x + ç fc 

qui implique que p divise q. D 

Enfin, il est nécessaire de pouvoir préciser les valeurs de h pour lesquelles la loi de 
Benford est satisfaite au moins lorsque les coefficients Cj sont rationnels. 

Théorème 6 Soit (a n ) n >i C M + \ {0} une suite qui satisfait la relation de récurrence 
d'ordre k 

0"n+k — Cfc-l a n+fc-l + Ck-2 a n+k-2 + • • • + C\a n+ i + CQa n 

avec des coefficients Ci > ; Cq > et {i < k : Cj > 0} ^ 0, et avec les valeurs initiales 
ai, . . . , Ofe strictement positives. Alors : 

(1) La suite (a n ) n >i satisfait la loi de Benford dans presque toute base b au sens suivant : 

lim — 1{2 < b < N : (a n ) ne suit pas la loi de Benford en base b}\ = 0. 

(2) Si les coefficients Q sont tous rationnels et si la racine réelle positive du polynôme 
caractéristique de la relation de récurrence n'est ni un entier ni un inverse d'entier, 
alors (a n ) n >i satisfait la loi de Benford dans toute base b > 2. 

Remarque. Dans le cas où seul Cq > 0, c'est-à-dire si / = 0, toute suite (a ra ) n >i qui 
satisfait la relation de récurrence (*) est de la forme : 

k-i 

^ = (^r-^«/* /fc Vn. 

1=0 

Par pédiodicité de e 2nlin / k ^ l a suite (a n ) est donc réunion de k sous-suites (a m+ k n )n>i, 
< m < k, qui sont toutes géométriques de raison cq. 

Ainsi, si Co 7^ 1, chaque sous-suite (a m+ kn)n>i satisfait la loi de Benford dans presque 
toute base. Si en revanche c = 1, chaque sous-suite est constante et ne satisfait pas la loi 
de Benford. 

Les preuves utilisent principalement la matrice compagnon du polynôme p(x) et elles 
reposent sur la théorie de Perron-Frobenius des matrices à coefficients positifs ou nuls qui 
sont irréductibles ou primitives, et dont nous rappelons les principaux résultats dans le 
paragraphe suivant. Les preuves des théorèmes 2 et 3 se trouvent dans le paragraphe 3 et 
la preuve du théorème 6 dans le dernier paragraphe. 

2 Matrices irréductibles ; matrices primitives 

Nous rappelons ci-dessous les définitions et les résultats principaux de la théorie de 
Perron-Frobenius à propos des matrices irréductibles et primitives. Nos références sont 
d'une part le chapitre 8 de la monographie de Cari D. Meyer [7] et d'autre part les notes 
de J E Rombaldi |g]. 



Soit C G Mfc(R) une matrice k x k k coefficients positifs ou nuls (on note C > 0). Nous 
considérerons ici les valeurs propres complexes de C, c'est-à-dire ses valeurs propres en tant 
qu'endomorphisme de C fc . L'ensemble des valeurs propres est le spectre de C et sera noté 
c(C). Nous rappelons que le rayon spectral de C est 

p{C) = max{|A| : A G a(C)}. 

On dit que C est réductible s'il existe une matrice de permutation P telle que PCP^ 1 soit 
de la forme : 

(/-il r-illi 
c» 

où C" et C" sont des matrices carrées de dimensions positives, et on dit qu'elle est irréduc- 
tible si elle n'est pas réductible. 

On démontre qu'une telle matrice C = (qj) est irréductible si et seulement si le graphe 
orienté G(C) associé est fortement connexe. Le graphe G(C) a pour ensemble de sommets 
les entiers 1, 2 . . . , k, et il y a une arête orientée de i vers j si et seulement si Cij > ; le 
graphe est fortement connexe s'il existe un chemin orienté de 1 vers 1 passant par tous les 
sommets. 

Si C est irréductible, le rayon spectral p = p(C) est strictement positif et c'est une valeur 
propre simple du polynôme caractéristique pc{x) = det(C — xi). Le sous-espace propre 
correspondant est de dimension (complexe) égale à 1, et il est engendré par un vecteur de 

Perron- Frobenius, c'est-à-dire un vecteur v = \ \ dont toutes les composantes t», sont 

\v k ) 

strictement positives. 

Notons h > 1 le nombre de valeurs propres A G o~(C) telles que |A| = p. Si h > 1, on 
l'appelle l'indice d'imprimitivité de C, et si h — 1, on dit que C est primitive. 

On démontre que si C est irréductible, alors elle est primitive si et seulement s'il existe 
un entier m > tel que C m > 0, c'est-à-dire les coefficients de C m sont tous strictement 
positifs (test de primitivité de Frobenius, cf [7], pp. 674 et 678). 

Si C est irréductible mais non primitive, si {Ai,..., A?,,} est l'ensemble des valeurs 
propres de C de module égal à p, alors 

{Ai, . . . , \ h } = {p, pu, pu 2 , ..., pu 11 ' 1 } 

oùu = e 2ni / h ([7], p. 676). 

Soit S = (sjj) G Mfe(M) une matrice à coefficients positifs ou nuls. On dit que S est 
stochastique (par rapport aux lignes) si, pour tout 1 < i < k, on a 



£ 



S i,j 



Si une telle matrice S est de plus irréductible, si h est son indice d'imprimitivité, alors elle 
admet toutes les racines /i-ième de l'unité comme valeurs propres. 

Grâce au test de primitivité de Frobenius, on observe que si A et B sont des matrices 
à coefficients positifs ou nuls irréductibles et s'il existe e > tel que A > eB (c'est-à-dire 
si aij > ebij pour tous i,j), alors A est primitive si B l'est. 



Considérons encore une matrice C > irréductible et notons v 



f v,\ 



\v k J 



> un 



vecteur de Perron- Frobenius de C ; la transposée C de C possède les mêmes propriétés, 



/w,\ 



donc elle admet un vecteur de Perron- Frobenius w 



> associé également au 



\w k J 
rayon spectral p. 

Si de plus C est primitive, alors la suite de matrices (4rC n )n>i converge vers la matrice 
G := vw T /w T v qui est la projection sur le noyau de C — pi parallèlement à l'image de 
C — pi Q7J, p. 674). On observe que G a tous ses coefficients strictement positifs. 



3 Preuves des théorèmes 2 et 3 
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Cfc-2 


1 

Cfc-l J 



Considérons d'abord un polynôme p(x) = x — Ck-\X — ■ ■ ■ — C\X — Cq comme dans 
la première partie du théorème 1, c'est-à-dire tel que Cj > pour tout 1 < j < k — 1 et 
cq > 0. On lui associe sa matrice compagnon 



C 



Cela signifie que le polynôme caractéristique pc(%) = det(C — xi) satisfait : pc(%) = 
(—ï) k p(x). Les deux polynômes ont par conséquent exactement les mêmes racines avec les 
mêmes multiplicités. Comme Cq > 0, C est irréductible puisque, dans le graphe associé, il y 
a une arête de 1 vers 2, de 2 vers 3, etc. de k — 1 vers k et de k vers 1, au moins. Ainsi, p{x) 
satisfait les conditions de la première partie du théorème 2 si et seulement si sa matrice 
compagnon est primitive. 

La partie (1) du théorème 2 est alors une conséquence immédiate de la proposition 
7 ci-dessous, qui est un cas particulier du théorème de la page 679 de [7J, mais nous en 
donnons une démonstration par souci d'être complet. 

Proposition 7 Soient cq, Ci, . . . , c^-i des nombres réels positifs ou nuls, Cç, > et soit 
I = {1 < j < k : Cj > 0}. Alors la matrice 



G 



( 1 
1 



y co ci 

est primitive si et seulement si pgcd(J U {k}) 



\ 





Cfc-2 Cfc_i J 

-- 1. 



Preuve. Comme nous l'avons observé ci-dessus, C est irréductible. Notons iV le cardinal de 



/. Soit S la matrice 



S 



( 






V n+î 



1 








1 





1 



o \ 





1 



■ N+l • • • N+l 

où dans la dernière ligne, pour tout j E I U {0}, Cj est remplacé par jA^- La matrice ainsi 
obtenue est donc stochastique. En particulier, p(S) — 1, et S est irréductible. De plus, par 
les rappels du paragraphe 2, C est primitive si et seulement si S l'est. On remplace alors 
C par S, qui est la matrice compagnon de 



q(x) 



x 



E 

je/u{o} 



x J 



N + l 



Supposons d'abord que d = pgcd(J U {k}) > 1. Écrivons j = a^d et k = ad avec a,j, a G N* 
pour tout j G /, et considérons le polynôme 



r(x) = x a — y 



je/u{o} 



X 3 

N+l 



Alors q(x) = r(x d ), r(l) = 0, et pour tout ( G C, ( d = 1, on a q(() = r(l) = 0. Cela 
démontre que C n'est pas primitive. 

Réciproquement, supposons que pgcd(J U {k}) = 1, et soit ( une racine de q(x) telle 
que \(\ = 1. Par le paragraphe 2, on sait que ( est alors une racine de l'unité. Ecrivons 
donc 



C 



^inm/d 



avec 1 < m < d et pgcd(m, rf) = 1. L'égalité q(() = donne 



(N + i) e 2inkm / d 



5> 

J67 



2iTrjm/d 



et en prenant les modules 



iV + 1 



V^ e 2iirjm/d _|_ ^ 



< 



E 



^iirjm/dl 



+ 1 = N + 1 



qui implique que d divise jm pour tout j G /, donc que e 2l7ï J rn / d = l ) puis que e 2mkm ^ d = 1, 
donc que d divise également km. Ces conditions impliquent que d divise pgcd(J U {k}) ; 
comme pgcd(m, d) = 1, on obtient que ( = 1, et S est primitive. D 

Pour démontrer la seconde affirmation du théorème 2, considérons un polynôme à coef- 
ficients réels p(x) = x k — cx k ~ 1 — E^~ c i x ^ ^ c l ue c > 2 et c > X^=o l c il + 1- On observe 



d'abord que pour tout z G C tel que \z\ > 1, on a 



fc-2 

E 



c.-z J 



< \z\ k -\c-l). 



Cela implique facilement que p(c— 1) < 0, que p(c+l) > et que p(z) ^ pour tout \z\ — 1. 
Ainsi, p s'annule en un nombre réel 1 < c — 1 < p < c + 1. Enfin, soit q(z) = z 



cz 



fc-i 



On a |ç(^)| > c — 1 pour tout \z\ = 1, et, pour ces mêmes valeurs de z, on a 

fe-2 



|p(z)-g(z)| 



£ 

3=0 



dz 3 



<c-l<\q(z)\. 



Puisque q admet k — 1 zéros (comptés avec multiplicités) dans le disque unité ouvert, il en 
est de même pour p par le théorème de Rouché. □ 

Considérons ensuite une suite (a n ) n >i comme dans la troisième partie du théorème 2 : 
a n > pour tout n, et elle satisfait la relation de récurrence 



l-n+k 



Ck-iûn+k-i + Cfe_ 2 a n+ fc-2 + • • • + cia n+ i + c a n (n > 1) 



où les Ci satisfont l'une des deux premières conditions du théorème 2. On note encore p 

. . — c de la 



x 



Ck-lX 



fe-1 



l'unique racine positive du polynôme caractéristique p{x) 
relation de récurrence. 

Pour démontrer la troisième partie du théorème 2, il suffit de vérifier que linin^oo ^ 
existe et est positive ([5], théorème 2.4). Le cas où les coefficients q satisfont la condition 
(1) est une conséquence de la proposition suivante. 

Proposition 8 Avec les hypothèses et les notations ci-dessus, si les Ci satisfont la condition 

(1) du théorème 2, on a : 

(i n 



lim 

n— >oo {?' 



existe et est positive. 



Preuve. Notons encore C la matrice compagnon du polynôme p(x) et posons pour tout 
entier n > 1 

\ 



Ar 



( a n 



\ O-n+k-1 J 

de sorte que 

CA„ = A n +\ 

pour tout n > 0. Par suite, on a A n = C n ~ x A\ pour tout n. Si v et w sont des vecteurs de 
Perron- Frobenius pour C et C T respectivement, on a vu que la suite de matrices (\C n ) 
converge vers la matrice G = vw T /w T v lorsque n — ¥ oo. Ainsi, 



P" 



-A n 



P P 



-C n ~ v A x -» 



p 



En particulier, ^, qui est la première composante de -^A n , converge vers la première 

composante de -G Ai, et comme Ai > et G > 0, toutes les composantes de GA 1 sont 
positives. □ 

Pour terminer la preuve du théorème 2, considérons une suite (a n ) n >i C M + \ {0} qui 
satisfait la relation de récurrence 



O-n+k — Ck-lCL-n+k-l + • • • + Cid n+ i + CoO r , 



avec Ck-i > 2, c^-i > J^ 7 =o \ c j\ + 1? e t telle que a n ne converge pas vers lorsque n — > oo. 
D'après la structure de telles suites (rappelée dans le paragraphe 1), si (j = p > 1, £ 2 , ■ ■ ■ , Cm 
sont les racines du polynôme caractéristique p(x) = x k — Ck-\X k ~ x — ... — C\X — Cq, de 
multiplicités respectives \x\ = 1, p 2 , • • • , Pm, il existe des coefficients a^i tels que 

a n = a 1A p n + Y^ Yl a M nl Q Vn > L 

j=2 1=0 

Puisque p > 1 et \Q\ < 1 pour tout j = 2, . . . ,m, a n est de la forme a n = a\^\p n + (3 n avec 
n — > lorsque n — >■ oo. Les hypothèses a n > pour tout n et a n 7^ impliquent que 
«i : i > 0, et on obtient immédiatement que linij^oo ^ = a^\. D 

Nous passons enfin à la preuve du théorème 3. 

Preuve du théorème 3. En n'écrivant que les coefficients non nuls dans l'expression du 
polynôme p(x), on a 

p( x ) = x fc - c fc _ fcl x fc " fci - ... - c k ^ ks x k ' ka - c 

avec 1 < &i < . . . < k s < k. Par le théorème de la page 679 de [7J, on obtient : 

h = pgcd(fc -ki,...,k-k s ,k)= pgcd(/ei, . . . , k s , k) = pgcd(7 U {k}). 

Ecrivons, pour 1 < i < s, ki — hk' { , et aussi k = hk! . Pour chaque m fixé, en remplaçant n 
par m + hn dans la relation de récurrence, on obtient 

a m+h(n+k') = c fc-fci G m+h(n+fc'-/4) + • • • + Ck-k s a m+h(n+k' -k' s ) + Co a m+/in- 

Cela signifie que la suite (a m+ h n )n>i satisfait une relation de récurrence d'ordre k' dont le 
polynôme caractéristique q(x) est 



/ \ k' k' — k' k 1 — k 1 

q(x) = x K - c k _ kl x k fc i - . . . - Ck-kX fcs 

avec pgcd(J(ç) U {k'}) = 1. On applique alors les conclusions des propositions précédentes 
à la suite {a m+hn ) n >x. D 

4 Critère d'irrationnalité de log^(p) 

Soit p(x) le polynôme p(x) = x k — Ck-\X k ~ l — ... — C\X — Cq avec Cj G M + et cç, > 0. On 
désigne encore par I l'ensemble des indices 1 < i < k tels que q > 0, et par p la racine 
positive de p(x). On suppose encore que 1^0. 

Le résultat suivant sera utilisé dans la preuve du théorème 6. 

Proposition 9 Soit b > 2 un entier. On suppose que les Ci G Q + et que : 

(1) pgcd(7 U {k}) = 1 ; 

(2) p n'est ni entier, ni inverse d'entier. 
Alors log 6 (p) est irrationnel. 



Preuve. Supposons que log 6 (p) soit rationnel, de sorte qu'il existe deux entiers p et q > 

premiers entre eux tels que p = b p ' q . Comme p ^ 1, on a p / 0. Désignons par m(x) le 

polynôme minimal (unitaire) de p. 

Si p est irrationnel, alors m(x) est de degré au moins 2 et p en est une racine simple. 

Mais p est également une racine de t(x) = x q — If, donc le polynôme m(x) divise t(x), et 

ainsi toutes les racines de m(x) sont des racines de t(x). Or, ces dernières sont toutes de 

module égal à p. Cela contredit (1) car m(x) divise également p(x), et ce dernier aurait 

plusieurs racines de module égal à p. 

Par suite, p est nécessairement rationnel ; il existe des entiers positifs a, (3 tels que 

p = § et pgcd(a, (3) = 1. On obtient a q = b p f3 q , et l'unicité de la décomposition en facteurs 

premiers implique que 

/3 = 1 si p est positif, c'est-à-dire si p > 1, et alors p = a est entier, 

a = 1 si p est négatif, c'est-à-dire si p < 1, et alors p = 1//3 est un inverse d'entier, 

ce qui contredit la seconde hypothèse. □ 

Remarque. La condition (2) est facile à utiliser puisqu'il est commode de localiser p dans 
M. + grâce à une étude succinte du polynôme p(x). 

Exemple 10 Soit m > 1 un entier fixé; pour tout entier k>2, soit 

Pk,m(x) = x k — mx k ~ 1 — mx k ~ 2 — ... — mx — m, 

qui généralise le polynôme caractéristique des suites de Fibonacci d'ordre k introduites dans 
l'exemple 4, et qui satisfait les conditions du théorème 2. Notons p^ m la racine positive de 

Pk,m(x). 

Elle est irrationnelle par la proposition 5, pk, m > 1 car p& !TO (l) = 1 — km < 0, et en 
fait, pk, m est un nombre de Pisot (cf [2J, [lj) : c'est un entier algébrique, et Pk,m{x) est son 
polynôme minimal car on a c^-i > ■ ■ ■ > Cq > 0, et toutes les racines ( ^ p^ t7n de Pk, m {x) 
satisfont |C| < Pk,m- 



On a m < pk :m < m + 1 car 



Vk,m\ m ) — \ m-m k 
m— 1 



k < 0, m 



< 0, m > 2, 



et Pfc >m (m + 1) = 1 pour tout m. 

Plus précisément, on va démontrer que p feîn < Pfc+i, m pour tout k assez grand et pour 
tout m, et que 

(m + l)k 

, 1 < Pk,m < m + 1. 

Cela démontrera que, pour tout m fixé, la suite (pk,m) est croissante et converge vers m + 1. 
La première affirmation provient des égalités : 

Pk+l,m(Pk,m) = Pk^m ~ m Pk,m ~ ■■■ ~ mp k , m ~m = p fc>m ■ p k ,m(pk,m) — W, = ~m 
et Pk+l,m{pk+l,m) = 0. 

Pour démontrer la seconde affirmation, nous introduisons le polynôme auxiliaire 
q(x) — (x — l)pk,m(x) = x h+1 — (m + l)x k + m. 
On vérifie sans peine que 

f(m+l)k\ m + 1 f(m + l)k\ k 

q \ k + i )-—kTï\ k + i ) +m - 
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Or, l'inégalité -f±± (^£0) + m < est équivalente à (m + l) fc+1 > m(k + 1)(1 + l/k) k 
qui est vraie dès que k est assez grand. Cela démontre la seconde affirmation. 

Remarques. (1) Les suites (p% m )t>i ne sont pas équidistribuées mod 1 car la distance 
entre {p e km : £ < n) et N tend vers 0. Mais qu'en est-il de (a • p f km )i>\ pour a irrationnel 
positif? (On sait que, pour presque tout nombre réel x > 1, la suite (x n ) n >i est équidistri- 
buée mod 1; cf [5J, chap.l, corollaire 4.2.) 
(2) Si Co, Ci, ... , Cfc-i G Z avec k > 2 sont tels que 

fc-2 

Cfc-i > y^|cj| + 1 

3=0 

et Co 7^ 0, le polynôme associé p(x) = x k — c^-xX l — ... — Cq est le polynôme minimal 
d'un nombre de Pisot p k ([lj, chap. 5.2), et cela donne une famille de suites qui satisfont 
une relation de récurrence à coefficients entiers non nécessairement positifs et qui satisfont 
également la loi de Benford d'après le théorème 2. 

Nous passons enfin à la preuve du théorème 6 : 

Preuve du théorème 6. (1) Pour N > 3 fixé, posons 

B N = {2 < b < N : log 6 (p) G Q}. 

Il suffit de démontrer que 



y/iVlog(iV) 
|5iv| " log(2) • 
C'est évident si B^ est vide ou s'il ne contient qu'un élément. S'il contient au moins deux 
éléments, soit &o — min{6 G -B/v}- Pour tout b G B^ tel que 6 > bo, en utilisant encore la 
décomposition en facteurs premiers, on vérifie qu' il existe un entier u < yN tel que 6 e ^ 
b soient des puissances entières de u. Par suite, 



B N C {2 < u p < N : p > 1, u < ^/N}. 



On obtient ainsi la majoration annoncée. 
(2) Soit h > 1 le nombre de racines positives distinctes du polynôme p(x). Les sous-suites 
qui ont des conditions initiales positives parmi {a m +hn)n>i, < m < h, ont un polynôme 
caractéristique dont p l l h est la racine positive. Si p n'est ni entier ni inverse d'entier, alors 
p l / h non p] us e t les conditions de la proposition 9 sont satisfaites. 

Enfin, si p est rationnel, en utilisant la décomposition en facteurs premiers, on vérifie 
comme ci-dessus que b ±p / q G N. Par suite, il existe des entiers c, d > 1 et u > 2 tels 
que p = u c et b = u d , mais alors p serait entier ou inverse d'un entier, contrairement à 
l'hypothèse. Cela démontre la deuxième partie du théorème 6. D 
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